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THE SEVEN DIMENSIONAL PERFECT DELAUNAY 
POLYTOPES AND DELAUNAY SIMPLICES 

MATHIEU DUTOUR SIKIRIC 


Abstract. For a lattice L of a sphere S{c,r) of center c and radius r is 
called empty if for any v ^ L we have Iju — c|| > r. Then the set 5(c, r) fl L 
is the vertex set of a Delaunay polytope P = conv(5(c, r) fl L). A Delaunay 
polytope is called perfect if any affine transformation 4> such that 4>{P) is a 
Delaunay polytope is necessarily an isometry of the space composed with an 
homothety. 

Perfect Delaunay polytopes are remarkable structure that exist only if n = 1 
or n > 6 and they have shown up recently in covering maxima studies. Here we 
give a general algorithm for their enumeration that relies on the Erdahl cone. 
We apply this algorithm in dimension 7 which allow us to find that there are 
only two perfect Delaunay pol3d:opes: 821 which is a Delaunay pol3d.ope in the 
root lattice E7 and the Erdahl Rybnikov polytope. 

We then use this classification in order to get the list of all types Delaunay 
simplices in dimension 7 and found 11 types. 


1. Introduction 

A lattice L is a set of the form L = Zti + • • • + Zz;„ C M" with (vi,... ,Vn) being 
independent. For such L a sphere S(c, r) of center c and radius r is called empty if 
for any v £ L we have ||ri — c|| > r. A polytope P is called a Delaunay polytope if 
it is full-dimensional and if the vertex-set of P is ^(c, r) Cl L with 5'(c, r) an empty 
sphere. A Delaunay polytope is called perfect if any affine transformation </> such 
that (j){P) is a Delaunay polytope is necessarily an isometry of the space composed 
with an homothety. 

In [25l [9] it was proved that for dimension n < 5 the only possible perfect 
Delaunay polytope is the interval [0,1]. Also in [9] it is proved that the Gosset’s 
polytopes 221 and 821 which are Delaunay polytopes of Eg and E7 are perfect. 
From the construction of infinite sequences in [m [221 ESI HI] we know that for any 
dimension n > 6 there exist perfect Delaunay polytopes. In |15| for any n > 6 we 
define a Delaunay polytope EDn of a lattice LD„. The lattice LDn is formed by 
lamination over the root lattice D„_i and we prove in [22] that EDn is the unique 
Delaunay polytope of maximum circumradius of LDn and compute its covering 
density. 

In [To] we proved that the EDq = 221 is the unique perfect Delaunay polytopes 
in dimension 6. This work uses a new approach in order to prove the following 
theorem: 

Theorem 1. The 7-dimensional perfect Delaunay polytopes are the Cosset polytope 
EDf = 821 and the Erdahl and Rybnikov polytope ERy [241128] . 
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Proof. See Section HI □ 

Perfect Delaunay polytope are of importance for the theory of Covering Maxima. 
A covering maximum is a lattice L such that its covering density is reduced if it 
is perturbed. In [22) it is proved that a lattice L is a covering maximum if and 
only if the Delaunay polytopes of maximum circumradius are perfect and eutactic 
(see [22] for the definition). This characterization echoes Voronoi’s theorem [44] 
for the characterization of lattices of maximum density in terms of perfection and 
eutacticity. In |22j we proved that LDn is one such covering maxima. Based on 
TheoremjTjand partial enumerations in dimension 8, 9 and 10 we state the following 
conjecture: 

Conjecture 1. For each n > 6, the lattice LDn defined in [m has maximal cov¬ 
ering density among all covering maxima. 

The Minkowski conjecture |341 p. 18] on the product of inhomogeneous forms 
has inspired a lot of research. Recently, it has been proved for n < 8 in mm 
Eol 129] by computational methods based on Korkinc-Zolotarev reduction theory. 
Other theoretical approaches have been attempted in [ST] [42] by Dynamical System 
Theory. In particular the following theorem is proved in jH] Corollary 1.3]: 

Theorem 2. If ConjectureUl holds for a dimension n >1 then Minkowski’s con¬ 
jecture holds for dimension n. 

As a consequence of the work of this paper we have that Minkowski’s conjecture 
is correct in dimension 7, thereby confirming |29j . 

We prove Theorem [1] by using the Erdahl cone which is defined as the set of 
polynomial functions / of degree at most 2 such that f{x) > 0 for x € IP. We 
already used this cone in |22j for the study of covering maxima. We have then 
to do a kind of dual description computation with the problem that the number 
of defining inequalities is infinite, we have no local polyhedrality result as in the 
perfect form case (see [41] for details) and we are interested in only a subset of the 
extreme rays (see Theorem [7] for the list of possible kinds of extreme rays of the 
Erdahl cone). 

In [To] we used a different approach, i.e. hypermetrics that allowed us to find all 
the 6-dimensional perfect Delaunay polytopes. But this approach relied on previous 
work [IHH] on 6-dimensional Delaunay simplices that we could not extend easily to 
dimension 7. Thus, it appears that the only way to classify the perfect 7-dimensional 
Delaunay is to use the Erdahl cone. Moreover, we are able to use this classification 
in order to get the classification of Delaunay simplices: 

Theorem 3. Up to arithmetic equivalence there are II types of 7-dimensional 
Delaunay simplices. The full list is given in Table [7] 

Proof. See Section [THl □ 

In contrast to perfect Delaunay polytopes, the lattices simplices of this list (ex¬ 
cept the trivial simplex) had not been discovered before. A similar study has been 
undertaken in [53] for the set of shortest vectors of lattices. In view of this work it 
seems reasonable to think that the classification of Delaunay simplices is possible in 
dimension 8. Equally importantly the classification of perfect Delaunay polytopes 
in dimension 8 could be done and a conjectural list of the 27 known possibilities is 
available in [16] . 


THE SEVEN DIMENSIONAL PERFECT DELAUNAY POLYTOPES AND DELAUNAY SIMPLICES 


Table 1. Representative of Delaunay simplices in dimension 7. 
ei,..., 67 is the standard basis of and cq = 0. vol(S') is n! times 
the Euclidean volume of S. |Stab(5')| is the size of the lattice 
automorphism group preserving S. “Nb interval” is the number of 
Delaunay polyhedra of the type {0,1} x Z® in which Si is contained 


i 

Representative Si 

vol(S'i) 

|Stab(5,)| 

Nb interval 

1 


eo, . . . , 66,67 

1 

40320 

127 

2 

Co, 

..., 66 ,( 1 , 1 , 1 , 1 , 1 , 1 , 2 ) 

2 

40320 

63 

3 

Co, 

...,e6,(0,0,l,l,l,l,2) 

2 

1440 

63 

4 

Co, 

...,66,(1,2, 2, 2,2, 2,3) 

3 

5040 

42 

5 

Co, 

...,66,(1,1,1,2, 2, 2,3) 

3 

1152 

42 

6 

Co, 

...,66,(0,!,1,2, 2, 2,3) 

3 

240 

41 

7 

Co, 

...,66,(1,1,1,1,1,3,4) 

4 

1440 

27 

8 

Co, 

...,66,(1,1,1,1,2, 2,4) 

4 

240 

31 

9 

Co, 

...,66,(1,1,1,2,2, 3,4) 

4 

144 

31 

10 

Co, 

...,66,(1,1,3, 3,3,4,5) 

5 

72 

24 

11 

Co, 

...,66,(1,1,1,1,2,3,5) 

5 

48 

24 


In Section [5] Delaunay polyhedra are considered, their basic structure and re¬ 
lation to the Erdahl cone are introduced here. The facial structure of the Erdahl 
cone is reviewed in Section [3l in particular not all extreme rays of the Erdahl cone 
are related to Delaunay polyhedra [25]. We also explain how the hypercube [0,1]" 
correspond to the cut polytope in the Erdahl cone. Section |4] is not used in later 
sections. In it we construct a retraction of the Erdahl cone on the faces defined 
by Delaunay polyhedra. In Section [5] we establish the link between the Erdahl 
cone and the classic T-type theory. In Section [6] we do the same for the hyper¬ 
metric cone. In Section [3 we give the connectivity and finiteness results on which 
our enumeration algorithm relies. Then we present in Section [8] our enumeration 
method, which is modelled on the Voronoi algorithm for perfect forms [35] and on 
the adjacency decomposition method |5]. In Section |9] we give the obtained results 
in the classification of 7-dimensional perfect Delaunay polyhedra. In Section [TOl we 
use this classification to classify the 7-dimensional types of Delaunay simplices. 

2. Delaunay polyhedra 

Denote by E 2 (n) the vector space of polynomials of degree at most 2 on R" and 
by AGL„(Z) the group of affine integral transformations on Z". The Erdahl cone 
is defined as 

Erdahl{n) = {/ g £’ 2 ( 71 ) such that f{x) > 0 for a; G Z"}. 

It is a convex cone of dimension {n + l)(n -I- 2)/2 on which the group AGL„(Z) 
acts. All defining inequalities f(x) > 0 are equivalent under AGL„(Z) and therefore 
Erdahl(n) is not polyhedral. 

We denote • the standard scalar product on M" defined hy x ■ y = x'^y. For a 
symmetric matrix A and x G R" we define A\x\ = x'^Ax. We write any / G E 2 {n) 
in the form 


f{x) = Gst(/) -H 2 Lin(/) • x + Quad(/) [x] 
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with Cst(/) € K., Lin(/) € M" and Quad(/) a n x n symmetric matrix. We define 
S'” to be the set of symmetric matrices, S”o the set of positive definite matrices, 
S”g the set of positive semidefinite matrices. We also define Erdahl^o{n) to be the 
set of / G Erdahl{n) with Quad(/) G S”g. If / G Erdahl{n) then Quad(/) G S^q- 
A sublattice of Z” is a subgroup of Z”. An affine sublattice is one of the form 
xq + L with Xq G Z” and L a sublattice of Z”. A lattice L C Z” is called saturated 
if (L(g)]R) nZ" = L. If Li and L 2 are two sublattices of Z” we write Li 0x^2 = 
if Ai nL 2 = {0} and Li UL 2 generates Z" over Z. In that case both Li and L 2 are 
saturated. 

Definition 1. Let us fix n > 1 and define: 

(i) A Delaunay polyhedron D to be a set of the form D = Pli{D) + L{D) C Z” 
with 

• conv(PL'(D)) a Delaunay polytope of an affine sublattice L' of , 

• L{D) a sublattice ofWP and 

• L'©zL(D) =Z”. 

(ii) A Delaunay simplex set S to be a Delaunay polyhedron with |S| = n + 1. 
(Hi) A repartitioning set R to be a Delaunay polyhedron with |S| = n + 2. 

The isotropy lattice L{D) is uniquely determined by D and its dimension is called 
the degeneracy rank of D denoted degrk(D). Note that D is the vertex set of a 
convex body only when L(D) = 0. Also Delaunay polyhedra is full-dimensional, 
i.e. the smallest affine saturated lattice containing D is Z” itself. 

The set Pl'{D) is included in L' and depends on L'. For any two choices Tj 
and L '2 there exist a bijective affine map (j) : L( ^ L '2 with (j){PLi^{D)) = P^{D). 
When we consider properties that do not depend on the integral representation, we 
drop the lattice and write P{D). 

For / G Erdahlin) we write 

Z{f) = {a; G Z” such that f{x) = 0} . 

In classical geometry of numbers the essential tool is the quadratic form Q instead 
of the quadratic function. The following establish a direct link between both: 

Definition 2. For a Delaunay simplex set S G IP and Q G S’^ there exists a 
unique function f G E 2 {n) such that: 

• f{x) = 0 for X G S 

• and Q = Quad(/). 

This function is denoted fs^Q and depends linearly on Q. 

The key reason for using Delaunay polyhedra is the following theorem: 

Theorem 4. (i) If D is a Delaunay polyhedron then there exist a function f G 
Erdahl{n) such that D = Z{f). 

(ii) If f G Erdahl{n) then either Z{f) is empty or there exist a k-dimensional 
saturated affine lattice L C Z” such that Z{f) is a Delaunay polyhedron of L. 

Proof, (i) Let us take a Delaunay polyhedron D = Pli{D) + L{D) with Pli[D) 
being a Delaunay polyhedron of a lattice L' with L' ©z L{D) = Z”. Let us denote 
by S(c, r) the sphere around Pl'{D) and write L' = Zui + • • • + Zu^. The function 

/' : Z” ^ R 

x = {xi,...,Xn) - c||^ - 
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belongs to Erdahl{k). More precisely, f'{x) = 0 if and only if ^ Pl'{D). 

For X S Z" we write x = xi+z with xi € L' and z G L{D) and write f{x) = /'(xi). 
It is easy to prove that / G Erdahl{n) and Z{f) = D. 

(ii) This is [25l Corollary 2.5]. □ 

We define the rational closure >q to be the set of positive semidefinite forms 
whose kernel is defined by rational equalities. 

Corollary 1. If f € Erdahl(n) is such that Z{f) is a Delaunay polyhedron then 
Quad(/) G 

Proof. Let us write D = Z{f) and take a lattice L' C IP with L' 0z L{D) = Z". 
We write any x G Z" as x = xi +z with Xi G L' and z G L{D). There is a quadratic 
function /i on L' such that /(x) = /i(xi) and Z{fi) = Pli{D). Thus Quad(/i) is 
positive definite and since L{D) is an integral lattice the matrix Quad(/) belongs 

to □ 

Given a set 1^ C Z" we will need to be able to test if it is a Delaunay polyhedron 
or not. Algorithm [T] does this iteratively for a finite point set by solving larger 
and larger linear programs until conclusion is reached. The algorithm can be easily 
adapted to the case of a point set of the form R + L with L a lattice and R finite. 
The corresponding algorithm for perfect form is given in [331 Algorithm 1]. 

If D is a n-dimensional Delaunay polyhedron, then we define 

Aut(D) = {(?!> G AGL„(Z) : (f{D) = D}. 

When using Algorithm [T] it is best to impose that the sought function / is invariant 
under Aut(D) since it simplifies the search and a Delaunay polyhedron admits an 
invariant function (see Corollary [2| below). 

Before stating our result on the description of Aut(D), we remind on the notion 
of semidirect product. Given a group G, we call G a semidirect product and write 
G = N X H N \s a, normal subgroup oi G, H a subgroup G = NH and N Cl H = 
{e}. 

Theorem 5. If D is a n-dimensional Delaunay polyhedron of degeneracy degree d 
then we have the isomorphism 

Aut(D) = ((Z'^)^+"-'^ XI GLd(Z)) X Aut(P(i:>)). 

Proof. Let us take a basis v = (yi)Ki<d of L{D). Any automorphism of D will 
send V to another basis of L{D) and this determines a component GLd(Z) of the 
automorphism group. Let us write L' for an affine sublattice of Z" such that 
L' 0z L{D) = IP. Such affine sublattice are determined by I 0 n — d vectors in 
L[D) and this determines the component (X<^)i+’^-d of the automorphism group. 
The last component comes from the fact that the automorphisms have to preserve 
the polytope conv(F(D)). □ 

We denote by Aff(D) the normal subgroup xi GLd(Z) given in the 

above theorem. 

Corollary 2. A Delaunay polyhedron D admits a function f G Erdahl(n) with 
Z{f) = D that is invariant under Aut(D). 

Proof. Let us write L' 0z L{D) = Z". A vector x G Z" is decomposed as x = 
xi + z with xi G L' and z G L{D). Any function / with Z(f) = D must have 
Quad(/)[z] = 0 and Lin(/) • z = 0. Therefore, there exist a function fi on L' such 
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Data: A finite set F C Z” of full affine rank 

Result: A quadratic function / G ErdahlyQ{n) such that Z{f) = 1/ if it 
exists and false otherwise 
Syert ^ V ■ 

Svect ^ 0- 

repeat 

V <— random element of Z" 

Syert Syert U 

until The set {evy for € Syert} has rank {n + l)(n + 2)/2; 
while no solution has been reached do 
Form the linear program 

minimize/gB^ln) Tr(Q,ua.d{f)) 

subject to f{v) = 0 for u S F 

f lv) > 1 for U e Syert - V 
Quad(/)[?;] > 1 for u G 

Syect 

if The linear program is infeasible then 
I return false 
end 

/o a rational optimal solution, 
if /o G Erdahl{n) and Z{fo) = V then 
I return /o 

end 

if Quad(/o) ^ S'>o then 

Find a vector u G Z" with Q{fo)[v] < 0 
Syect t Syect U 

end 

if Quad(/o) G and /o ^ Erdahl{n) then 
Find a vector u G Z" with /o(u) < 0 
Syert t Syert U{v} 

end 

if /o G Erdahl{n) and Z{fo) ^ V then 
Find a vector v G Z{fo) — V. 

Syert ^ Syert U {u} 

end 

end 

Algorithm 1: Testing Delaunay realizability of a finite set of points 


that f{x) = fi{xi). So, / is invariant under Aff(D). The group acting on /i is 
Aut{PL'{D)), which is finite. The function 

fiixi) = /(^(a^i)) 

u^Aut{Pj^f {D)) 

is Aut(PL'(D)) invariant. Thus we get a function f'{x) = f'i{xi) invariant under 
Aut(D). □ 

If D and D' are two Delaunay polyhedra such that D C D' then we define the 
stabilizer group 

Stab(D, D') = {(j}G AGL„(Z) : = D and = D'} 

= Aut(D) n Aut(D'). 
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We have the following results 


Theorem 6. Let D and D' be two Delaunay polyhedra satisfying D <Z D'. Then: 

(i) we have L{D) C L{D') and AfF(I?) C Aff(iD'). 

(ii) There exist a finite group Gi C Aut(P(-D)) such that 


Stab(£),£)') = Afr(D) x Gi 


in particular Aff (D) is a finite index subgroup of Stab(D, D') 

Proof. We have L{D) C L[D'). Let us take a Z-basis {ei,..., Cd} of L{D) and com¬ 
plement it by adding {cd+i,..., e^'} to a basis of L{D'). We can find /i,..., fn-d' 
so that the Ci and fj form a basis of Z". The group Aff(£)') is generated by transla¬ 
tions along /i,..., fn-d' and GLd'{Z). The group generated by translations along 
{cd+i, ■ ■ ■ ,ed'} and GLd(Z) directly embeds into GLd'CZ) and this determine the 
group inclusion. So (i) holds. 

By (i) we have the inclusion AS{D) C Stab(D,Z?'). Let us choose a lattice 
L' with L' ©z L{D) = Z”. Then if m S Stab{D, D') C Aut(D) we can find an 
unique element n € Aff(£)) such that un~^ stabilizes L'. Thus un~^ belongs to 
Aut{PL'{D)). The image defines the group Gi. the finite index property follows. 


□ 


3. Facial structure of the Erdahl cone 


The standard scalar product on S'” is {A,B) = Tr{AB). We equip Erdahl{n) 
with the inner product 

U.g) = Cst(/) Cst(g) + 2Lin(/) • Lin(g) + {Quad(/), Quad( 5 )) 


and for each x € Z” we define the evaluation function evx by evx{y) = {1 + x • y)^ 
such that 


(f,evx) = f{x). 


A convex cone C is defined as a set invariant under addition and multiplication 
by positive scalars. C is called full-dimensional if the only vector space containing 
it is R™. C is called pointed if no linear subspace of positive dimension is contained 
in it. Let C be a full-dimensional pointed convex polyhedral cone in R™. Given 
/ G (R™)*, the inequality f{x) > 0 is said to be valid for C if it holds for all x € C. 
A face of C is a pointed polyhedral cone {x G C : f{x) = 0}, where f{x) > 0 is a 
valid inequality. 

A face of dimension 1 is called an extreme ray of C; a face of dimension m — 1 
is called a facet of C. The set of faces of C forms a partially ordered set under 
inclusion. We write F <l G if F C G and dimF = dimG — 1. Two extreme rays 
of C are said to be adjacent if they generate a two-dimensional face of C. Two 
facets of C are said to be adjacent if their intersection has dimension m — 2. Any 
(m — 2)-dimensional face of C is called a ridge and it is the intersection of exactly 
two facets of C. 

By the Farkas-Minkowski-Weyl Theorem (see e.g. SOI Gorollary 7.1a]) a convex 
cone C is polyhedral if and only it is defined either by a finite set of generators 
{wi,..., uat} C R™ or by a finite set of linear functionals {/i,..., /m} G (R’”)*: 
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Every minimal set of generators {z;i,..., vn' } defining a polyhedral cone C has 
the property 

= {e : e extreme ray of C}. 

Every minimal set of linear functionals {fi,, /m'} defining C has the property 
that {Fi,... ,Fm'} with Fi = {x G C : fi{x) = 0} is the set of facets of C. The 
problem of transforming a minimal set of generators into a minimal set of linear 
functionals (or vice versa) is called the dual description problem. 

In our work, we have to deal with Delaunay polyhedra with an infinite number 
of vertices and we cannot apply the Farkas-Minkowski-Weyl theorem to them nor 
of course existing dual-description software mm- 

Definition 3. Let D C IP be a Delaunay polyhedron 

(i) We define the vector space 

Space{D) = {/ G E 2 {n) such that f{x) = 0 for x G D}. 

(ii) The dimension of Space{D) is called the perfection rank rankperf(Z?) and 
D is perfect z/rankperf(Z?) = 1. 

Proposition 1. Let D C IP be a Delaunay polyhedron 

(i) Space{D) n Erdahl{n) is a face of Erdahl{n) of dimension dim Space{D). 

(ii) If D is perfect then Space{D) n Erdahlin) is an extreme ray of Erdahlin). 

Proof. Let p = dimSpace{D) and gi, ..., Pp he a basis of Space{D). Since D is 
a Delaunay polyhedron there exist a function / G Erdahl{n) such that D = Z{f). 
Eor each 1 < i < p there exist Ai > 0 such that Xif + gi G Erdahl{n) and so (i) 
follows, (ii) follows directly from (i). □ 

For a perfect Delaunay polyhedron D we denote fo a generator of the extreme 
ray Space{D) n Erdahl(n). 

Theorem 7. f |25l Theorem 2.1]^ The generators of extreme rays of Erdahl(n) are: 

(1) The constant function / = 1 

(2) The functions of the form fa.pix) = (oiXi . a„a;„ -1-/3)^ with (ai, ..., a„) 
not colinear to an integral vector. 

(3) The functions of the form fu with D a perfect Delaunay polyhedron. 

This theorem indicates that the structure of the extreme rays of Erdahl(n) is 
more complicated than for a polytope. Since we are interested only in the third 
class of extreme rays, some reduction will be necessary and it turns out that we 
can work out everything with Delaunay polyhedra. 

In this paper we will work with both spaces of functions in E 2 (n) and with point 
sets of Delaunay polyhedra: 

Definition 4. Given a Delaunay polyhedron D we define 

(i) The cone of admissible functions is defined as 

Erdahlsupp{D) = {g G ^ 2 ( 71 ) : g{x) > 0 for all x G D} . 

(ii) The cone of evaluation functions is defined as 

Erdahl*^pp{D) = E XxCVx with Ax > 0 
I xGD 
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Theorem 8. Let D be a Delaunay polyhedron of perfection rank r. 

(i) Erdahlsupp{D) is the product of a pointed convex cone Cd with Space{D). 
The dual of Cd is Erdahl*^^^{D). 

(ii) Any Delaunay polyhedron D' <Z D of perfection rank r + 1 gives a facet of 
Erdahl*^pp(D). 

(Hi) If L(D) = 0 then facets of Erdahl*^pp{D) correspond to Delaunay polyhedra 
D' C D of perfection rank r + 1. 


Proof, (i) If / and — / both belong to ErdahCuppiD) then f{x)=0toixGD and 
so / G Space{D). So, Erdahlsupp{D) is the sum of Space{D) and a closed convex 
cone. The duality result follows from [31 Part IV.5]) for closed full-dimensional 
convex cones. 

(ii) Let / G Erdahl{n) such that Z{f) = D'. Thus we have f{x) =0 on D' and 
f{x) > Q on D — D' and so / defines a facet of Erdahl*^pp{D). 

(iii) If L{D) = 0 then for any facet of Erdahl*.^pp{D) we can find a set D' C D 

and a function / G E 2 {n) with f[x) = 0 for a; G and f{x) > 0 for x G H — D'. 
There exist a function g G Erdahl{n) such that D = Z{g). Then we can find A > 0 
such that f + Xg G Erdahl{n). Then we have D' = Z{f -\- {\-\- 1 ) 5 ) and so D' is a 
Delaunay polyhedron. □ 

Theorem |25l Theorem 2.1] shows that if L{D) ^ 0 there are other facets of 
Erdahll.^pp{D) than the ones from Delaunay polyhedra. 


Proposition 2. If D is a n-dimensional perfect Delaunay polyhedron with degen¬ 
eracy degree d then P{D) has at least 




points. 


Proof. The isotropy lattice L{D) has dimension d and we can choose a complement 
lattice L' such that L' ©z L{D) = Z”. If / is a quadratic function with Z{f) = 
D then / is determined by its restriction to L'. Hence / belongs to a vector 
space of dimension (”” 2 '*'^) gives the minimal number of determining 

inequalities. □ 

A surprising relation has been found between the Erdahl cone of the hypercube 
{0,1}" and the cut polytope, which is classic polytope of Combinatorial Optimiza¬ 
tion [12]: Write = {1,..., n}; if S' C A^, then the cut metric Ss on N is defined 
as follows: 

if |{*j}ns| = l, 

I 0 otherwise. 

We have 63 = Sn-s and the cut polytope CUTP„ is defined as the convex hull of 
the cut metrics 83 . The cone defined by the cut polytope is defined as 


CCUTP„ = ^ As(l, 83 ) with As > 0 

y 5c{l,...,n} 

The facets of the cone CCUTPjj are in one-to-one correspondence with the facets 
of the polytope CUTP„. 


Theorem 9. The polyhedral cone Arda/iZ*„pp({0,1}”) is linearly equivalent to 
CCUTP„+i. 
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Proof. The hypercube [0,1]" is defined as the convex hull of 2” vectors v = {vi,... ,Vn 
with Vi G {0,1}. For every such vector the evaluation function is 

evy(x) = 1 + 2 (v ■ x) + (v ■ xY- 
Thus we can associate to ev^ the vector 

(1; (u2)l<2<nj (UiUj )l<j<j<7i). 

Since vf = Vi for Ui S {0,1} this vector family is linearly equivalent to 

(1; ) l<i<j <n ) ^ ^ {0, 1} ■ 

To the same v we can associate the vector v = {0,vi,... ,Vn) and the set S = 
{i G {0,...,n}, s.t. zJi = 1} C (0, The cut metric 5$ on {0, ...,n} is 

characterized by 

j))o<i<j<n 

and since Ss{i,j) = (vi —VjY the family of pairs (l,i5s) is linearly equivalent to 
the family of evaluation map evy D 

It is interesting to note that the symmetry group of the hypercube [0,1]" is of 
size 2"n! but that the symmetry group of the cut polytope CUTP„+i is of size 
2"(n+ 1)! for n 3 [HI 

4. The Delaunay polyhedra retract 

For / G Erdahl{n), we define VectZ{f) to be the vector space spanned by 
difference of elements of Z{f). We define Vf = VectZ{f) + Ker Quad(/). For 
a given / G Erdahl(n) we call proper pair a pair (g^h) G E 2 {ri)‘^ such that g G 
Erdahl(n), h{x) > 0 for a: G M" and f = g + h. 

Lemma 1. Let f G Erdahl{n). 

(i) For a proper pair {g, h), Z{f) C Z(g) and Ker Quad(/i) C Vf. 

(a) There exist a proper pair {g, h) with Ker Quad(/i) = Vf 

Proof. If a: G Z{f) then one sees that necessarily h{x) = g{x) = 0. So, h{x) = 0 
for X G conv(2'(/)) which implies VectZ{f) C Ker Quad(/i). Also, it is clear that 
for any vector v G Ker Quad(/) we have Quad(g)[u] = Quad(/i)[z;] = 0. Hence, (i) 
holds. 

Let us denote by ZKer Quad(/) the smallest subspace of M" having an integral 
basis containing Ker Quad(/). By the Decomposition Lemma 3.1 in [25] there exist 
a g G Erdahl(n) with Ker Quad( 5 ) = ZKer Quad(/) and a positive semidefinite 
form Qi such that 

f{x) = g{x) + Qi[x]. 

Let us denote by V, respectively W, an integral supplement of Ker Quad(/) in 
ZKer Quad(/), respectively ZKer Quad(/) in Z". Denote by ^i, ..., fim some 
affine functions on Z" such that (/)i(ZKer Quad(/)) = 0 and 

jx G IF s.t. = • • • = (j)m{x) = 0} = Vect Z{f) n W. 

Then for e > 0 small enough the function gi — is still in Erdahl{n). 

So, one gets that the pair (/ — h, h) with h(x) = Qi[x\ + is proper 

and (ii) is true. □ 

Let us call IF an integral supplement of Vf. Denote by Quad(/)|vT the qua¬ 
dratic form Quad(/) restricted to IF. A proper pair {g, h) is called extremal if 
det Quad(/i)|w is maximal among all proper pairs. Lemma |TJ(i) implies that the 
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notion of being extremal is independent of the chosen subspace W, while Lemma 
[T](ii) implies that there is at least one form of non-zero determinant. 

Theorem 10. Let f € Erdahl{n). 

(i) An extremal proper pair has Z{g) being a Delaunay polyhedron. 

(ii) There exist a unique extremal proper pair [g, h). 


Proof. Let us take an integral supplement W as above and suppose to avoid triv¬ 
ialities that IT 7 ^ 0. So, by restricting / to W, we can assume that Quad(h) is 
positive definite. 

f(x) = g(x) + h(x) 


with h(x) > 0 for a; G 


‘. This condition on h is equivalent to 
Cst(h) Lin(/i)^ 


Ah = 


Lin(h) Quad(h) 


being positive definite. Hence we consider the following semidefinite programming 
problem: find the Ah £ <S'>g ^ maximizing det Quad(h) and satisfying for all a; G Z” 


f{x) > Ah[{l,x)\ = h{x). 


We also write g = f — h. 

Suppose that Z{g) is not a Delaunay polyhedron and that Quad(/i) is positive 
definite. Then Z{g) does not generates R" as an affine space and so there exists an 
affine function </> such that 4>{Z{g)) = 0. Then there exist a > 0 such that the pair 
(/ — h', h') with h' = h + aff' is still proper. Since det Quad(/i') > det Quad(/i) 
the pair is not extremal and (i) holds. 

Let us take N = n{n+l)/2 points Vi G IP such that the family {(1, Vi){l, 
is of full rank. The inequalities g{vi) > > 0 implies that all coefficients 

of A are bounded. Thus the problem is actually to minimize the convex function 
h I—>■ — logdet Quad(h) over a compact convex set, hence existence follows. 

Since — log det is a strictly convex function, we know that if we have two optimal 
solutions hi and h 2 then Quad(/ii) = Quad(/i 2 ). Let us denote by Di = Z{gi) 
and D 2 = Z{g 2 ) the corresponding Delaunay polyhedra. The function hmid = 
{hi + h2)/2 is also an optimal solution of the problem. We have Z{gmid) = L)i nD 2 . 
The set Di r\D 2 is necessarily a Delaunay polyhedron since otherwise we could still 
increase the determinant by the above construction and this would contradict the 
optimality. But if Z{f) is a Delaunay polyhedron then the terms Cst(/) and Lin(/) 
are determined by Quad(/). So, one gets hi = /12 and the uniqueness is proved 
on the restriction to W. But Lemma [DJi) implies that once Quad(h) is known on 
W then it is known on Z”. By the condition h > 0 the linear part is known as 
well. □ 

Note that in the above determinant maximization problem a finite set of in¬ 
equalities suffices to determine the optimal solution. This follows from the fact 
that since we are maximizing the determinant we can assume that the lowest eigen¬ 
value of Quad(h) is bounded away from 0, i.e. that there exist c > 0 such that 
Quad(/i) > cin. 

For / G Erdahl(n), we write proj{f) = g and proj'{f) = h with {g,h) the 
unique extremal pair associated to /. From the unicity of extremal pairs we also 
get that proj and proj' commute with the action of AGL„(Z). 


Conjecture 2. The function proj is continuous. 
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Let US define Erdahldp{n) to be the set of / G Erdahl{n) such that Z{f) is a De¬ 
launay polyhedron. The above conjecture if true implies that the set Erdahldp{n) 
is simply connected and this could be of interest for topological applications. How¬ 
ever, we were not able to prove the conjecture and instead we prove in later sections 
connectedness results which are sufficient for our purposes. 

5. Relation with L-types theory 

In this section we reframe classical L-type theory from [3^ (see also [H] for a 
modern account) in term of Erdahl cone and state several key lemmas. 

Definition 5. Let Q G >q. The Delaunay polyhedra tessellation 'DPT{Q) de¬ 
fined by Q is the set of Delaunay polyhedra D such that there exist a f G Erdahl (n) 
with 

• Zif) = D 

• and Quad(/) = Q. 

If Q is positive definite then the Delaunay polyhedra tessellation is the classical 
Delaunay polytope tessellation, i.e. all Delaunay polyhedra occurring are actually 
vertex sets of Delaunay polytopes. The number of translation classes of Delaunay 
polyhedra is always finite. These Delaunay polyhedra tessellations were considered 
in [20l Section 2.2]. Efficient algorithm for the enumeration of Delaunay polytope 
tessellations are given in [21] . 

From this one can define the L-type which are parameter spaces of Delaunay 
polytope tessellations: 

Definition 6. Let us take a Delaunay polyhedra tessellation T. Then the L-type 
LTfT) is defined as the closure of the set of quadratic forms Q such thatDPT{Q) = 
T■ It is well known (see [41) and |45) for proofs) that L-types are polyhedral cones. 

A L-type is called primitive if it is of maximal dimension, this is equivalent to 
say that all its Delaunay polyhedra are Delaunay simplex sets. 

The set of all L-types for all possible Delaunay tessellations defines a tessellation 
of the cone 

Given two Delaunay polyhedra tessellation E and E' we say that E' is a re¬ 
finement of E if every Delaunay polyhedron of E' is included in a single Delaunay 
polyhedron of E. E' is a simplicial refinement if all its Delaunay polyhedra are 
Delaunay simplex sets. 

Proposition 3. Any Delaunay polyhedra tessellation E admits at least one sim¬ 
plicial refinement. 

Proof. Let us denote by L{E) the space L{D) of the Delaunay polyhedra D occur¬ 
ring in the tessellation and by Q G >q the form realizing it. Let us take a lattice 
L' such that L' 0z L{E) = IE. We write x G IE a.?, x = x\ P z with x\ G L' and 
z G L{E). D is a Delaunay polyhedron for the quadratic function / G Erdahlfn). 
Necessarily / is of the form f{x) = fi{xi) with fi a quadratic function on L'. Let 
us denote {Di)i^i the Delaunay polyhedra occurring in the tessellation. Let us take 
a basis wi, ..., Wm of L{E) and define linear forms fi on IE such that fifwj) = Sij 
and 4>i{L') = 0. The quadratic form 

m 

Q'[x\ = Q[x\ + ^ 

2 = 1 
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is positive definite and the Delaunay polyhedra tessellation VT corresponding to 
Q' is formed by the Delaunay polyhedra 

m 

Di + Ufc + llwfe with 1 < i < r and 

k=l 

In particular 72 is a refinement of 'T. 

Since the i-type domain form a tiling of >q the form Q' belongs to at least 
one primitive L-type LT(72)- This L-type defines a Delaunay polyhedra tessellation 
by simplices which is a refinement of 72 and so of T. □ 

Let us take a primitive L-type T- Any facet A of T is determined by a pair 
of Delaunay simplex sets Si and S 2 in the Delaunay tessellation that determine 
a repartitioning set. We say that two facet dehning repartitioning sets are in the 
same class if they define the same facet F of T. If i? is a repartitioning set then 
conv(i?) admits exactly two triangulations (see [HJ Section 4.3.2]). One says that 
two primitive L-types are adjacent if their intersection is a codimension I face in 
the cone S'"o- When we move from one L-type to another L-type, the Delaunay 
tessellation is changed and this is done combinatorially by the repartitioning sets. 
That is some Delaunay simplex sets are merged into repartitioning sets and the 
triangulation is changed to the other triangulation, thus yielding another L-type. 

Given a Delaunay polyhedron D a Delaunay polyhedra tessellation T is called 
D-proper if D is the union of the Delaunay polyhedra D' contained in D. We have 
following lemma: 

Lemma 2. Let D be a Delaunay polyhedron. The graph formed by the primitive 
L-types whose corresponding Delaunay polyhedra tessellations is primitive and D- 
proper is connected. 

Proof. Let us consider a function fjy £ Erdahl{n) such that Z{f]j) = D. We can 
consider the triangulations induced by positive definite quadratic forms on D itself. 
This set is connected by the theory of regular triangulations (see [7] for an account). 

Any triangulation Fpart on D induced by a positive definite quadratic form Q can 
be extended to a triangulation T of Z": It suffices to replace Q by Q-\- X Quad(/D) 
for A sufficiently large. The reason is that Quad(/£)) will not change the Delaunay 
triangulation for Delaunay simplex sets contained in D. 

Now given a primitive L-type LT whose Delaunay polyhedra tessellation F is 
D-proper, we denote by S its set of Delaunay simplex sets included in D. We 
consider the following cone C{S): 

C{S) = {g £ s.t. fs,Q{x) > 0 for S' £ 5 and a: £ Z” - 5} . 

This cone is convex and is an union of primitive L-types. Thus this set of L-types 
is connected. The connectedness follows by combining above results. □ 

Delaunay polytopes restricted to the set D. When doing so, by the connectivity 

Theorem 11. Let us take D a Delaunay polyhedron and two Delaunay simplex 
sets S and S' in D. Then there exist a sequence of Delaunay simplex sets {S = 
So, Si,..., Sm = S'} with Si C D for 0 < i < m such that Si U S^+i is a reparti¬ 
tioning set for 0 < i < m — 1. 

Proof. Let us take fo £ Erdahl{n) a function such that Z{fo) = D. Take fs, fs' 
the corresponding functions for S and S'. Denote by T the Delaunay polyhedra 
tessellation dehned by Quad(/£)), which obviously has D as one of its component. 
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Table 2. Possible volume of lattice Delaunay simplices. See [5] 
for the proof for n < 6 and Section [10] for the proof for n = 7. 


n 

1 

2 

3 

4 

5 

6 

7 

PossVol{n) 

{1} 

{1} 

{1} 

{1} 

{1,2} 

{1,2,3} 

{1,2,3,4,5} 


When we perturb Quad(/£)), we are changing the Delaunay tessellation. How¬ 
ever, if we take e > 0 small enough, we can ensure that the Delaunay polyhedra 
tessellations X>PT(Quad(/£) + e/s)) and X>PT(Quad(/D + ef'g)) are D-proper. By 
applying Proposition |3| we can find simplicial refinement of those two tessellations 
which we name TR and TR' and are both D-proper. We call LT and LT' the 
corresponding primitive L-types. 

By Lemma [2] there exist a path between LT and LT' that uses only D-proper 
L-types. By following this path, we can change S into another Delaunay simplex 
set S 2 in TR' ■ 

Denote by /i, ..., fr the facets of LT'. Every such facet corresponds to a 
family of repartitioning sets. We say that two Delaunay simplex sets included in D 
are adjacent if their union is a repartitioning set which gives a facet of LT'. The 
Delaunay polyhedron D is a coarsening obtained by merging all simplices, so the 
above defined graph is connected. This means that we can find a path from S 2 to 
S'. □ 


6 . Relation with hypermetric theory 


We define the volume vol(S') of a Delaunay simplex set S to be n!PoZ(conv(5')) 
with Vol the Euclidean volume. This rescaled volume is an integer and satisfies 
vol(S') < n\. The possible rescaled volumes PossVol(n) are given in Table |6| for 
n < 7 and a super-exponential lower bound on max PossVol{n) is proven in [39) . 
The best known upper bound [121 Proposition 14.2.4] is 


( 1 ) 


max PossVol(n) < n! 


2 " 



Definition 7. Let us take two n-dimensional Delaunay polyhedron D, D' with 
D <Z D'. We can define the generalized hypermetric cone 

Hyp(D, D') = {/ € £’ 2 ( 74 ) s.t f{x) = 0 if x € D and f{x) >0 if x G D'} . 

We have the inclusion Hyp(D, D') C Erdahlsupp{D') and Hyp(Il, D') is a priori 
defined by an infinity of inequalities. 

As a direct application we can express the L-type domains as intersection of 
generalized hypermetric cones: 

Proposition 4. Let T be a Delaunay polyhedra tessellation. Then we have 

LT{T) = n^er QuadHyp(D,Z"). 

Proposition 5. The polyhedral cone Hyp(Zl, D') is polyhedral. 

Proof. Let us take a Delaunay simplex set S = {uo,...,Un} C D which exists 
by Proposition [3l If we prove the polyhedrality of Hyp(S', D') then Hyp(D,D') is 
polyhedral as well since it is obtained from Hyp(S', D') by adding equalities f{x) = 0 
iov X G D — S. 
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Suppose that a v G D' defines a relevant inequality. Then there exist a function 
/ such that /(x) = 0 for x G S' U {w} and /(x) > 0 for x G D' — S U {x}. 
Since D' is a Delaunay polyhedron, there exist a function g such that g{x) = 0 
for X € D' and g{x) > 0 for x G Z” — D'. Then we can find A > 0 such that 
f{x)+Xg{x) > 0 for X G Z"—SU{x}. As a consequence the polytope conv(SU{x}) is 
a Delaunay polytope. This implies that for any i G {0, ... ,n} the Delaunay simplex 
set Sv^i = {iijUo,... ... ,Vn} has vol(S„,i) < n\ (see proof of Theorem 

14.2.1 in [H]). Hence the coefficients of v are bounded by a bound depending only 
on S and this proves that Hyp(S, D') is polyhedral. □ 

For a Delaunay simplex set S of volume 1, the cone Hyp(S, Z”) is called the hy¬ 
permetric cone and is studied in [12) . For other simplices, they are called Baranovski 
cone in [41j . The facets of the Baranovski cones are determined up to dimension 
6 in [38]. There is a correspondence between facets of Hyp(S', D) and repartition¬ 
ing sets P with S C P C D. That is the inequality /(x) > 0 defines a facet of 
Hyp(S', D) if and only if S' U {x} is a repartitioning set. 


7. Connectivity results 
For a given Delaunay polyhedron D, let us write 



D' C D s.t. D' Delaunay polyhedron 
rankperf(D') = r and dimL(Zl') < d 


If A and B, are sets of Delaunay polyhedra then the graph Gr{A, B) is the graph 
on A with two Delaunay polyhedra Di,D 2 G A adjacent if and only if Di nD 2 G B. 

Theorem 12. If D is a Delaunay polyhedron of perfection rank r and degeneracy 
degree d then Cj.+i^d{D) is decomposed into a finite number of orbits under Aut(Il). 

Proof. Without loss of generality, we can write D = P-l-Z'^. Let us write D' C D as 
D' = P' -\- L' with k = dim Lh By applying an element of Aut(Il), we can assume 
that L' = Z^. So, without loss of generality, we can assume that L' = 0. Let us 
take a Delaunay simplex set S in D'; its volume is bounded by max PossVol{n). 
Again by using Aut(I?) we can find a constant C" such that the absolute value of 
the coordinates of S' in P -|- Z'^ are bounded by C'. The polyhedrality of the cones 
Hyp(S, D) implies the finiteness. □ 

Lemma 3. If C is a polyhedral cone, F a face ofC and e, e' are two extreme rays 
which are not contained in F, then e and e! are connected by a path which does not 
intersect F. 

Proof. By taking the intersection C H H with H a suitable hyperplane, we can 
transform C into a polytope P and e, e' into vertices of P. We can find an affine 
function (j) such that (fix) > 0 is a valid inequality on P and (fix) = 0 defines the 
face F n H of P. By maximizing the function (f over P and using the simplex 
algorithm (see [lOl ESj), we can find paths piv,Vopt)^ pivfvopt) from v, v' to an 
optimal vertex Vopt such that f is monotone on both paths. Since fiv) > 0 and 
fiv') > 0, such paths avoid the face F and put together gives the required path. □ 

Theorem 13. If D is a Delaunay polyhedron of perfection rank r and degeneracy 
degree d > 1 then 


Gr (Cr+l,d(D), Cr+ 2 ,d-liD)) 


is connected. 
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Proof. Let us take two Delaunay polyhedra Da and in D of perfection rank 
r + 1 . Let us take two Delaunay simplex sets S, S' contained in Da, By 
using Theorem IIII we can find a chain of simplices (S'i)o<i<m with Si C D, Sq = S 
and Sm = S'. Denote by Ri = Si U Si+i the repartitioning set. Write e_i, Cm 
to be the extreme rays in Hyp(S'o,D), Hyp(5'm,D) corresponding to Da, D'^. For 
each Delaunay simplex set Si we consider the cone Hyp(S'i,L)). The extreme rays 
correspond to Delaunay polyhedra of rank r +1. If a Delaunay polyhedron D' G D 
has degeneracy degree d then necessarily L{D') = L{D). We define the restricted 
trace function to be 

</>(/) =T’r(|Quad(/)|^(£,)) . 

A function / G Erdahlsupp{D) has Z{f) Ci D of degeneracy degree d if and only if 
</>(/) = 0. The hyperplane (j){f) = 0 determine a face Fi of the cone Hyp(S'i,D). 
The intersection is 


Hyp(5'i,i:)) n Hyp(5'i+i,i:)) = Ryp{Ri,D). 

Thus we can find a ray in Hyp(i?i,D), which is not contained in Fi. Since 
Hyp(5'i,D) is polyhedral, by Lemma [3] there exists a path from e^-i to Ci in 
Hyp(5'i,D) that avoids the face Fi. So, by putting all the paths together, we 
got the required connectivity result. □ 

Lemma 4 . If Di and D^ are two Delaunay polyhedra of perfection rank r and r + 2 
with D3 C Di then there exist exactly two Delaunay polyhedra D2.1 and D2^2 with 

D3 C D2,i C Di 

with rankperf(Z32,i) = r + 1. 

Proof. Since D3 is a Delaunay polyhedron there exist a Delaunay simplex set S C 
D3. The Delaunay polyhedra D3, Di correspond to faces F3, Fi of dimension r + 2 , 
r in the cone Hyp(S', Z"). It is well known from polytope theory |461 Theorem 
2.7.(iii)] that there are exactly two faces ^ 2 , 1 , ^ 2,2 containing Fi and contained in 
F 3 . Those gives the corresponding Delaunay polyhedron. □ 

By using this theorem, we are able to compute inductively the Delaunay poly¬ 
hedra in Z”. The property with the degeneracy degree ensures that we are able to 
effectively reduce the complexity of the computation at each step and thus we are 
reduced in the end to computation with Delaunay polyhedra of degeneracy 0, i.e. 
polytopes for which polytopal methods exist. 

8 . Algorithms 

In [2 a general survey of methods for computing dual description of highly 
symmetric polytopes with many facets are presented. Among the method presented 
there, we want to adapt the Recursive Adjacency Decomposition Method to our 
situation, i.e. to a case with an infinite group and an infinity of defining inequalities. 

8 . 1 . Computing Aut(D). In this subsection we explain the techniques needed to 
compute Aut(ZI), compute Stab^i'(Aut(ZI)) and split orbits. In the decomposition 
of Theorem[5]the only component that is not clear is Aut(ZIi), i.e. the computation 
of the automorphism group of a Delaunay polytope. For that purpose the methods 
of [21j can be used. The one that we are using is the method of isometry groups. 
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Definition 8. Let D be a Delaunay polyhedron. Take a lattice L' with L'(Bi.L{D) = 
Z”. Take a basis wi,... ,Wr of L'. Denote by vi,... ,Vm the expression of the 
vertices of Pl'{D) in the basis (wi). 

We define the matrix Q by 



From then we define the distance function : Z” x Z" i—>■ K 6?/ 

• foix.x') = 

• with <f>{x) = (1, Ml, . . . , Ur) if 

X = uiwi + • • • + UrWr + 2 and z € L{D). 

The construction of the matrix Q and its inverse above is relatively standard. 
We used it first in |43j and further work on this are done in 011]. 

The interest of this construction is that it allows to compute automorphism 
groups. 

Theorem 14. Let D be a Delaunay polyhedron. It holds: 

(i) If u £ Aut(Il) then we have fD{u{x),u{y)) = foix^y) for x,y G D. 

(ii) If L' is a sublattice such that L' ©z L{D) = Z” and u is a permutation of 
Pl'{D) such that fD{u{x),u{y)) = foix^y) for x,y £ D then u is induced by an 
affine rational transformation of L' 

Proof, (i) By its construction if 2; G L{D) then we have foix + z,y) = fD{x,y). 
Thus if M G Aff(£)) then u preserve fo- On the other hand if u G Aut{PL'{D)) then 
we can see by summation that u preserve fo- The proof is available for example in 

010 . 

(ii) The reverse implication is also available from [4310. □ 

Let us denote by AutQ(P(Zl)) the group of rational transformations preserving 
P{D). By Theorem [m (ii) we have Aut(P(£))) = AGLr.(Z) n AutQ(P(Il)) with r 
the dimension of L'. The computation of AutQ{P{D)) is done efficiently by using 
known partition backtracking software such as |36| : see 00 for more details. In 
the cases considered in this paper the number of vertices is quite small and this 
computation is very easy. 

A Delaunay polytope is called generating if difference between its vertices gener¬ 
ate Z”. If a Delaunay polytope is non-generating, then it is actually a Delau¬ 
nay polytopes for more than one lattice. If P{D) is generating then we have 
Aut(P(D)) = Autq{P{D)) and we are done. Otherwise, we can apply some of 
the strategies listed in 0 Section 3.1]. Here, the situation is particularly simple 
and the simplest strategy of iterating over the group elements and keeping the in¬ 
tegral ones works very well. Also note that the above methods with only slight 
modifications work for testing equivalence of Delaunay polyhedra. 

8.2. Computing stabilizers. We now give methods for computing stabilizers of 
Delaunay polyhedra, more precisely of the transformations preserving two polyhe¬ 
dra D <Z D' which occurs in our computations. 

Let us select a lattice L' such that L' ®z L{D) = Z". Denote by Gi the group 
occurring in Theorem 0 

Let us define the following function on P^i {D) by 

fD.D'ix,y) = {fD{x,y),fD'{x,y)). 
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By Theorem 1141 the elements of the group Gi must preserve the function Jd.D'- 
Thus we may use the partition backtrack algorithm of Section l8T] to get the group 
Aut(/£) £i/) of permutations preserving /d.d'- 

Then we obtain the group Gi by keeping only the elements that are in AGL„(Z) 
and preserve D' . This is possible since the group Aut{fD,D') is finite and of mod¬ 
erate size in most cases. All the algorithms above have equivalents for testing 
equivalence and of course what has been done for pairs D C D' of Delaunay poly- 
hedra can be extended to triples D C D' C D". 

8.3. Splitting orbits. Suppose that we have an orbit Gx of an element x under a 
group G. For a subgroup H C G we wish to decompose Gx into orbits Hxi. Such 
an orbit splitting decomposition 

Gx = U^iHxi with Xi = giX 
is equivalent to a double coset decomposition 

G = U™iiJg, StabG(x). 

In the case of interest to us we have G = Aut(D), H = Aut(D) n Aut(D') for 
D, D' Delaunay polyhedra with D G D' and x a Delaunay polyhedron included in 
D. Since a priori Aut(D) is infinite we cannot apply standard tools from computer 
algebra software such as GAP [28]. By the finiteness result Theorem [101(11) we can 
find a coset decomposition 

G = 

However, it is not a double coset decomposition, i.e. we can have Hgi ^ Hgj but 
still have HgiX = HgjX. Therefore, we need to eliminate duplicate in order to do 
the orbit splitting. 


8.4. The flipping algorithm. Suppose that Di C D 2 C D 3 are Delaunay poly¬ 
hedra having 


rankperf Di = 1-1- rankperf D 2 = 2 + rankperf D 3 . 

By Lemma |4| we know that there exist a unique Delaunay polyhedron D 2 having 
Di C D 2 C D 3 and D 2 ^ 

We can find functions /i S Erdahl{n) such that Z{fi) = Di. We can also assume 
that Gd = R/i © Goi+i for * = 1,2. We need to find /2 G Erdahl{n) such that 
Z{f'2) = D' 2 . 

If L^D^) = 0 then D 3 is a polytope, i.e. it has a finite number of vertices and 
the algorithm is called the gift wrapping procedure 0133115]. If L{D 3 ) ^ 0 we have 
to modify the algorithm in order to take care of the fact that we have an infinity 
of vertices by writing an iterative algorithm. This is quite similar to the flipping in 
the Voronoi algorithm m- 

The non-negativity test for /2 on Z" is done by solving a closest vector problem. 
The non-negativity test on D3 is done by decomposing it into {vi + L{D3)} U ■ • • U 
{vm + L{D3)}. The non-negativity is tested by m closest vector problems. The 
final operation on is done to ensure that /^{x) > 0 if x ^ -^(/s)- 
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Data: Delaunay polyhedra Di, D2, D3 with Di = Z{fi), fi € Erdahl{n), 
-Di C £>2 C D3 and 

rankperf Di — 2 = rankperf D2 — 1 = rankperf D3. 
Result: € Erdahl{n), Delaunay polyhedra = Z{f2), D2 ^ D2 and 

DiC D'2C Da 

V^0 

repeat 

V ■<— random element of Z{fa) 

V^VU{v} 

C ^ {fi{w),f2{w)) for w e V 
until C has rank 2; 

repeat 

{(ofi, /?i), (0:2, /32)} t— generators of extreme rays of C. 
h ^ {ai/i +/5 i/2, 02/1 + ^2/2} -R/2. 
if there is a v € Z{fa) with /^(v) < 0 then 
Ve-VU{u}. 

jC e- (/i(u’),/2(w)) for w; e V 

end 

until /2 > 0 on Z{fa); 

repeat 

I /2 ^ /2 + /a 

until /2 > 0 onV^; 

/2 ^ /2 + h 

Algorithm 2: Flipping algorithm 


8 . 5 . The recursive adjacency decomposition method. Given a Delaunay 
polyhedron D of perfection rank r the algorithm of this section will give the orbits 
of Delaunay polyhedra D' C D oi perfection rank r + 1. If degrk(D) = 0 then 
the computation of the orbits of Delaunay polyhedra can be achieved by Theorem 
11 ( 111 ). 

If degrk(D) > 0 we have to proceed differently. By Theorem [13] we can limit 
ourselves to Delaunay polyhedra with degrk(D) < degrk(D) — 1. The algorithm 
takes one initial Delaunay polyhedron of perfection rank r + 1 and computes the 
adjacent Delaunay polyhedron of perfection rank r + 1. If an obtained Delaunay 
polyhedron is not equivalent to an existing one then we insert it into the list. We 
iterate until all orbits have been treated. The computation of the adjacent Delaunay 
polyhedra adjacent to a Delaunay polyhedron D' requires the computation of orbits 
the Delaunay polyhedra contained in D' . Thus we have a recursive call to the 
algorithm. Fortunately the degeneracy degree diminish by at least 1 so there is no 
infinite recursion. 

The mapping from Ei to J -2 is done using the orbit splitting procedure. Finding 
D" from D 2 , D' and D is done using the flipping procedure. 

There is a degree of choice in the initial Delaunay polyhedron Dinit- The stan¬ 
dard choice is if D = P{D) + L{D) with L{D) = Zui -I- • • • -I- ZU(jegrk(D) to take 
Dinit a Delaunay polyhedron of the form 

Dinit = P{D) + {0, Ui} -f 1a]2 -f • • • -f ZU(iegrk(D)- 
The schematic of the algorithm is shown in Algorithm jSj 
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Data: Delaunay polyhedron D of perfection rank r 

Result: Set F = {Di,, Dm} of Delaunay polyhedron of perfection rank 
r + 1 inequivalent under Ant D 
if degrk(D) = 0 then 

U ^ orbits of facets of Erdahl*^pp{D). 

T •<— orbits of Delaunay polyhedron from Theorem |S1 (iii). 
else 

T •<— {Dinit} with Dinit € ErdahKji), rankperf Z[Dinit) = r + 1 and 
(legvk{Dinit) = degrk(D) - 1. 

while there is a D' G T with degrk(D') < degrk(D) — 1 do 
F^FU{D'}. 

T^T\{D'}. 

F\ ^ orbits of Delaunay polyhedra of perfection rank r + 2 in D' 
under Aut D'. 

Fi •(— orbits of Delaunay polyhedra of perfection rank r + 2 in D' 
under Stab(D',Z)). 

for Di € Fi do 

find Delaunay polyhedron D" C D of perfection rank r + 1 with 
D2 =D' n D". 

if D" is not equivalent under Aut D to an element of F UT then 
if degrk(D") = degrk(D) then 
I F^FUif} 
else 

I T^TU{f} 

end 

end 

end 

end 

end 

Algorithm 3: Enumeration of inequivalent sub Delaunay polyhedra 


9. Perfect Delaunay polytopes in dimension 7 

In the enumeration of inhomogeneous perfect form in dimension 7 we need to 
describe the Delaunay polytopes that will occur. The list of perfect Delaunay 
polyhedra in dimension 7 is thus 

(1) {0,1} X Z6. 

(2) 221 X ^ with 221 the Schlafli polytope 

(3) 2 > 2 i the Cosset polytope [T^ . 

(4) ERj the polytopes discovered by Erdahl and Rybnikov PH US]. 

The geometry of the Schlafli and Cosset polytope are described in more details in 

[Hill]. 

An affine basis of a n-dimensional Delaunay polytope D is a family of n + 1 
vertices vq, ... ^ Vn such that for any vertex v oi D there exist Ai G Z such that v = 
Sr=o The perfect Delaunay polytopes of dimension 7 have an affine basis but 
it is possible that in higher dimension there are perfect Delaunay polytopes without 
affine basis. It is known that in dimension at least 12, there are Delaunay polytopes 
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with no affine basis m- Also the perfect Delaunay polytopes of dimension 7 
are generating. Note that in [18] we found some non-generating perfect Delaunay 
polytopes for n > 13. We have rankperf({0,1}") = n (see for example [HI I12j l. 

In terms of computation, the overwhelming majority of the time is spent com¬ 
puting the rank 2 faces of {0,1} x Z®. By the recursive approach chosen, the method 
requires the computation of the facets of Erdahl*^pp{{0, 1}^) and so by Theorem 
M of the facets of CUTPg. We actually computed the list of orbits of facets of 
CUTPg (and some other graph cut polytopes) in [TT]. In dimension 8 the partial 
enumeration algorithm of m found 27 perfect Delaunay polytopes and it is likely 
that the list is complete. But to prove its completeness by using the method of this 
work would require the determination of all facets of CUTPg and this is very hard 
[6]. In [13| a partial enumeration of perfect Delaunay polytopes was done with only 
Delaunay polyhedra with L{D) = 0 being considered. The two perfect Delaunay 
polytopes of dimension 7 were determined in this work and our enumeration proves 
that the list is complete. 

The implementation is available from m and uses the GAP computer algebra 
system [28] . 

10. Classification of Delaunay simplices in dimension 7 

Formula © gives 187 as an upper bound on the volume of Delaunay simplex sets. 
With this upper bound we can devise an algorithm for enumeration of Delaunay 
simplex sets, which will unfortunately prove inefficient: 

Remark 1. Suppose we have a list of types of Delaunay simplex sets in dimension 
n — 1. If S = {uoi ■ ■ • ,Vn-i,Vn} « Delaunay simplex set of dimension n then 
{uoi ■ ■ •; Vn-i} is a n — 1 dimensional Delaunay simplex sets of the lattice vi — vq , 

■ ■ ■ ; ^ n —1 ^0 ' 

If we have a n — 1 dimensional Delaunay simplex set vq = 0, Vi € Z"“^ for 1 < 
i < n — 1 and is of volume v. We write the n-dimensional simplex as Vq = (vg, 0), 
v[ = (ui, 0), ..., v'^_^ = (u„_i, 0) and v'^ = (xi,... ,x„_i,x„) with x„ > 0. The 
volume of the simplex defined by [v'f) is vxn- 

For a fixed x„ the number of possibilities for (xi,... ,x„_i) is x”“^/u. 

The Delaunay simplices of dimension 6 were classified in [2] and so we could 
apply the algorithm of Remark |T] Unfortunately the number of possibilities to be 
applied is very large, of the order of 187®, on which we have to apply Algorithm |T] 

Therefore, we need a different method: 

Lemma 5. Let S be a Delaunay simplex set that is not contained in any perfect 
Delaunay polyhedron different from {0,1} x Z"“^. The possibilities are: 

(1) For n < 4 the Delaunay simplices of volume 1. 

(2) For n = 5 the Delaunay simplex of volume 1 or 2. 

(3) For n = 6 or 1 there are no possibilities. 

Proof. Let us take a Delaunay simplex set S = {uoi ■ ■ • i Pn}- We can assume that 
Vq is located at the origin by using translation if necessary. For each 1 < i < n let 
us define £i the linear form on M" such that £i(vi) = 1 and £i{vj) = 0 ior i i. 
Any Delaunay polyhedron D isomorphic to {0,1} x Z"“^ and such that S C D 
corresponds to a linear form £ on R" such that £{vi) £ (0,1} and 

D = (x € Z” I £{x) = 0 or 1}. 
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The linear form t is then called admissible and the corresponding quadratic function 
is qe{x) = £(x)(£(x) - 1). 

Let us denote by S C {1,..., n} the set of points i such that £{x) = 1. Clearly, 
one can write 

ies 

A function £ is admissible if and only if £ is integral valued on Z”. If it is not 
integral valued then there exists a u € Z” such that 0 < £{v) < 1 which implies 
that qi{v) < 0 which is not allowed. If it is integral valued then D = Z{qe) is 
equivalent to {0,1} x . In the following, for a set S' C {1,..., n} we write 

£s = '^£i and vs = 

ies ies 

Let us define 

S = {S C {I,..., n} s.t. £s is integral valued} . 

Let us denote by ZS the Z-span of the elements vs for S G S. This defines a 
lattice £ of Z”. Since S is contained only in Delaunay polyhedra isomorphic to 
{0,1} X Z”“^ the set of the function q^^ is full-dimensional. This implies that 
|S| > n{n + I)/2 and that the lattice £ is actually full dimensional. Denote by h 
its index. 

The set of the function is is also full-dimensional in (Z")*. Its index is 

h/vo\{S) > 1. 

We are interested in the point sets of the form {0, l}"n£ with L an affine subspace 
of Z". By direct enumeration we obtain the full list of 3363 orbits of such points for 
n = 7. By selecting the point sets whose cone of functions q^^ is full-dimensional 
we get an upper bound of 3 on the index h and so an upper bound of 3 on the 
possible volumes of such simplices. 

With volume at most 3 we can apply the algorithm implied by Remark [TJ Each 
facet of such a Delaunay simplex is also a Delaunay simplex of one dimension lower. 
Therefore, we can use previous enumeration result to get a list of 796 possible 
candidates of 7-dimensional Delaunay simplices. We then use Algorithm [T] for 
checking which ones of them are indeed Delaunay simplices. This gives 6 cases (the 
ones of Table [T] of volume at most 3). Each one of them is also contained in a 
Delaunay polytope ER'i and so there is no such Delaunay simplices in dimension 
7. 

Dimension n < 6 follows from known results. □ 

Proof of Theorem [S] If S' is a Delaunay simplex set then Hyp(S, Z^) is a 
full-dimensional polyhedral cone, i.e. defined by a finite number of inequalities and 
having a finite number of extreme rays. Any such extreme ray corresponds to a 
perfect Delaunay polyhedron D. We have |S| = 8 and S defines a face of the cone 
Erdahl*^pp{D). 

By Lemma [3 S has to be contained in a Delaunay polyhedron of type 32i, ER-j 
or 221 X Zj. 

The perfect Delaunay polyhedron ERj has 35 vertices which matches the lower 
bound given by Proposition [2] As a consequence any 8-element subset of ER'j 
defines a face of Erdahl*^pp{ER'j). The automorphism group of ER^ has size 1440 
and by using it one can get easily the 9434 orbits of 8-element subsets of ER^. 
Actually all 11 types of simplices occur this way. 


THE SEVEN DIMENSIONAL PERFECT DELAUNAY POLYTOPES AND DELAUNAY SIMPLICES 


The Gosset polytope 821 has 56 vertices and the automorphism group is equal 
to the Weyl group of the root lattice E7. We found 521 orbits of 8-element sets in 
821, 474 of them correspond to faces of Erdahl*^pp{32i)■ 

For the perfect Delaunay polyhedron 221 x Z we have to proceed differently since 
the number of points to be considered is infinite. We have to enumerate the possible 
8-point subsets of 221 x Z of volume at most 187 (Formula IH)) up to the action of 
Aut(22i X Z). The 8 points are expressed in the form u, = {wi, hi) with wt £ 221 and 
hi £ Z. The set of points (?i’i)i<i<8 must define a 6-dimensional affine space. Thus 
7 of them, say (wi)i<i<7, must be sufficient to define a 6-dimensional Delaunay 
simplex set Ssch- An exhaustive enumeration on the 27 vertices of 221 gives 81 
types up to isomorphism. The volume vol{Ssch) can be 1, 2 or 8. If the volume 
is 1, then we can use an element of AfF(22i x Z) and obtain hi = 0 for 1 < i < 7. 
For higher volumes, the situation is more complicate but by using AfF(22i x Z) and 
linear algebra we we can reduce to yo\{Ssch) possibilities, i.e. 2 or 8. For the last 
point (ug, hs) we have 27 possibilities for ug and a finite number for hg due to the 
upper bound of 187. We then apply Algorithm [T] to test realizability of the finite 
list of possible cases. This gives us the 11 possible simplices. □ 
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